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We consider the dynamics of Bose-Einstein condensates in a corral-like potential. Compared to 
the electronic quantum corrals, the atomic quantum corrals have the advantage of allowing direct 
and convenient observation of the wave dynamics. Our numerical study shows that this advantage 
not only allows people to explore the rich dynamical structures in the density distribution but also 
makes the corrals useful in many other aspects. In particular, the corrals for atoms can be arranged 
into a stadium shape for the experimental visualization of quantum chaos, which has been elusive 
with the electronic quantum corrals. 
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I. INTRODUCTION 

Quantum corrals were first demonstrated by arrang- 
ing iron adatoms into a ring on a copper surface [l|, Q . 
Several years later, an optical analogy to the electronic 
quantum corrals was proposed theoretically and realized 
in experiment by a skillful arrangement of nanoscale pil- 
lars 0, However, in both of the quantum corrals only 
static wave properties were studied experimentally. Since 
the density image was obtained by sequential scanning 
over a period of time, the experimental study of wave 
dynamics in these corrals is difficult if not impossible. 

In this paper, we study the dynamics of atomic quan- 
tum corrals for Bose-Einstein condensates (BECs). Since 
the density distribution of a BEC can be imaged snapshot 
by snapshot with a charge coupled device (CCD), the 
atomic quantum corrals offer a great advantage over their 
electronic and optical counterparts, namely, the possible 
experimental study of the wave dynamics inside the cor- 
rals. Our numerical simulation shows that rich dynamical 
structures of a condensate can arise in the quantum cor- 
rals due to reflection and interference. Moreover, with 
the great deal of experimental control over the conden- 
sates [HI, Q , one is possible to study novel quantum be- 
haviors, which are unimaginable for electronic or optical 
quantum corrals. For example, the dynamic evolution of 
a quantized vortex confined inside quantum corrals can 
be studied. 

Of particular importance, the atomic quantum corrals 
proposed here can be used as a laboratory to study quan- 
tum chaos. Even since the first creation of quantum cor- 
rals, it has been pursued to build stadium-shaped quan- 
tum corrals to visualize experimentally the "scar" states, 
a signature of quantum chaos [7]. The stadium-shaped 
quantum corrals were built; however, the visualization of 
quantum chaos has remained elusive because the quan- 
tum corrals are too "leaky" for electrons [t*] . Our atomic 




a blue-detuned 
laser beam 



an opaque sheet 
^ with a ring of holes 



a lens 



a condensate in a 
corral-like potential 



FIG. 1: Scheme to build quantum corrals for a Bose-Einstein 
condensate (see text). 



quantum corrals are very flexible and can be made to 
minimize the "leakage" so that be used to explore quan- 
tum chaos experimentally. As an example, we demon- 
strate with numerical simulation that with atomic quan- 
tum corrals of stadium shape, one should be able to visu- 
alize experimentally the quantum chaotic behavior pre- 
dicted in Ref. [8]. We also show how the interference 
between two BECs can be destroyed by chaos. 

The paper is organized as follows. In Sec. II, we con- 
sider the wave packet dynamics for a BEC in a corral- 
like potential. In particular, the vortex evolution in a 
corral-like potential is studied. In Sec. Ill, we study the 
dynamic quantum chaos for a BEC in stadium-shaped 
quantum corrals. A brief summary and discussion is 
given in the last section. 



II. WAVE PACKET DYNAMICS IN 
CORRAL-LIKE POTENTIAL 



* Electronic address: 
"I" Electronic address: 



xionghongwei@ wipm . ac . cn] 
bwu@ aphy. iphy. ac . cn | 



We consider here a scheme to create quantum corrals 
for a BEC illustrated in Fig. 1. An opaque sheet with a 
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ring of holes is placed between a lens and a blue-detuned 
laser beam to create a corral-like repulsive potential for 
the condensate. The lens is to focus and adjust the sizes 
of the "corrals" . The holes in the opaque sheet are iden- 
tical with a radius around or larger than 50/im, which is 
large enough so that the diffraction of the laser beam can 
be ignored. This setup is similar to the one used to create 
a rotating quasi-2D optical lattice by a mask with a set 
of holes [9] and a random potential ^ for BECs. One 
may also create this kind of quantum corrals by spatial 
light modulator or microlens array [llj. 

To simulate the dynamics of the cigar-shaped BEG in 
the corrals, we can integrate out the axial degree of free- 
dom. In this case, the BEC dynamics is described by a 
two dimensional Gross-Pitaevskii (GP) equation [6|]: 

.d^ _ _ (d^ ^ 

dt \ dx^ dy'^ J 

Voc^^VM^^g2Dm^^. (1) 

where ^ is a normalized wave function. The above equa- 
tion has been made dimensionless with the length unit 
Lq and the energy unit Eq = Ji^ /2mLQ. The time unit is 
then To = H/Eq. Vqc is the corral-like potential while 
Vht = muj\Ll {x^ + y'^) /2Eq with being the har- 
monic frequency in x and y directions. The dimensionless 
coupling constant g2D = 2\^Nh'^a/mlzEoLQ, where a 
is the 5— wave scattering length and Iz = ^JhjmijJz with 
(jjz being the trapping frequency in the z direction. We 
numerically solve this GP equation to illustrate the wave 
dynamics inside the quantum corrals. In the numerical 
calculations, we consider 2 x 10^ Na atoms in the conden- 
sate and use Lq = 5 /im, Tq = 1.8ms, and cc;^ = 4 x 27r 
Hz. 

The ith Gaussian potential Vi (x, y^ z) created by a fo- 
cused laser beam propagating along z direction has the 
following form 

Vi{x,y,z) ^ -^^^^-\^--.f+(y-y,f\l-\ (2) 

TTCr^ \Z) 

where P is the total power of a laser beam, a {z) = 
(Jo a/1 + I ^\ with the Rayleigh length zr = 2i\g^ / A. In 
the following calculations, ctq = 5 /im. For a focused laser 
beam with wave length A = 500 nm, we have zr — 314 
/im, which is much larger than the length (about 10 /im 
for the parameters in this paper) of the cigar-shaped con- 
densate. In this situation, along the z direction, the con- 
densate feels almost the same corral-like potential in the 
x — y plane. Thus, the two dimensional GP equation can 
be used very well to study the dynamics of a condensate 
in the corral-like potential. 

To show clearly the fundamental properties of the 
atomic quantum corrals, we calculate numerically the 
evolution of the condensate confined in circle-shaped 
quantum corrals, which can be described as a series of 

Gaussian potentials Kc = £ 7e-((^-^^-)'+(^-^^-)')/^' 
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FIG. 2: Density distribution of a condensate in a corral- like 
potential at different times: (a) t = 0, (b) t = 1, (c) t = 2, 
(d) t = 3, (e) t = 4, (f) t = 5. The coordinates are in unit of 
Lq. In Fig. a, the corral-like potential is also shown. 

with {xj^yj\ distributed uniformly along a circle of ra- 
dius R. With the parameters 7 = 20, M = 20, = 10, 
(7 = 1, and cj^ = 70 X 27r, we first get the ground-state 
wave function by using the widely used ITP (imaginary 
time propagation) method. Then the evolution of the 
condensate is solved numerically from the GP equation 
([1]), after suddenly switching off the harmonic trap. In 
this situation, the condensate will evolve under the con- 
finement of the corral-like potential. In Fig. 2a, both 
the initial density distribution |<I>|^ and the corral-like 
potential are shown. In Figs. 2b-2f, the evolution of |$| 
is given, after suddenly switching off the harmonic trap. 
The rich and colorful structures in |$|^ originate from 
the following two physical mechanisms: 

(i) After switching off the harmonic trap, and with the 
expansion of the condensate, the condensate will be re- 
fiected by the corral-like potential. The expanding and 
refiected condensates will overlap, and lead to clear inter- 
ference patterns. In Figs. 2c-2f, a series of ring-shaped 
interference fringes are clearly shown. 

(ii) In Figs. 2c and 2d, the density distribution shows 
sunfiower-like structure. This sunfiower-like structure is 
due to the discrete characteristic of the corral-like poten- 
tial. For the expanding condensate, the corral-like poten- 
tial can be regarded as M scattering sources, arranged 
along a circle. This discrete characteristic gives impor- 
tant modulation on the interference fringes. Because of 
this, the ring-shaped interference fringe is composed of 
a series of small peaks. In Figs. 2c and 2d, the number 
of small peaks distributed along an interference fringe is 
found to be exactly 20. By varying M, we have verified 
numerically that the number of small peaks distributed 
along a ring-shaped interference fringe is always equal to 
M, the number of corrals. 

The coherent interaction between atoms and external 
field (such as laser and magnetic field) can provide us im- 



FIG. 3: Evolution of the density distribution in the presence 
of both a corral-hke potential and a harmonic potential. The 
dynamic evolution is due to the sudden decreasing of the har- 
monic frequency, (a) t = 0, (b) t — 7r/5, (c) t — 27r/5, (d) 
t = 7r/2, (e) t = 37r/5, (f) t = Ttt/IO, (g) 47r/5, (h) tt. Here 
the coordinates are in unit of Lq. Fig. i gives the time evo- 
lution of the standard deviation 5 of the density distribution. 
In Fig. a, the corral-like potential and harmonic potential 
are also shown. The height of these two potentials has been 
adjusted in Fig. a. 



portant opportunity to study various dynamic processes 
[5| . We consider here the evolution of the condensate af- 
ter the sudden decreasing of the harmonic frequency to 
cj_L = IT.Stt. The evolution of the density distribution is 
shown in Fig. 3. The density distribution in Figs. 3d-3f 
shows typical flower-like structure in the kirigami. Be- 
cause of the confinement of the harmonic trap, we see 
a quasi-periodic oscillation behavior. In Fig. 3i, the 
standard deviation 5 of the density distribution is given, 
which shows further this quasi-periodic behavior. With- 
out the corral-like potential and interatomic interactions, 
the time period would be exactly tt/ (cj^^o)- From Fig. 
3i, the time period is slightly smaller than tt/ (a;±To) be- 
cause of the corral-like potential and interatomic inter- 
actions. 

The corral-like potential plays a key role in the inter- 
ference fringes of the density distribution. Without the 
corral-like potential, the decreasing of the harmonic fre- 
quency will only lead to a width oscillation of the conden- 
sate [6]. We have verified this by numerical calculations, 
which show that there are no interference fringes for this 
case. 

We now consider the evolution of a quantized vortex in 
the circle-shaped quantum corrals. In our calculations, 
we use 7 = 20, M = 20, = 10, a = 1, andcj^ = 70x27r. 
The initial vortex state in the presence of both the corral- 
like and harmonic trapping potentials is obtained numer- 
ically from the ITP method with a trial wave function, 

which has the general form f [\/x^ 2) ei'" The 
signs '+' and '— ' in e^'^ denote the rotational directions 
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FIG. 4: Density distributions of a quantized vortex in a corral- 
like potential at different times. The left column is for a 
vortex rotating counterclockwise at (a) t = 2.5, (b) t = 3.8, 
(c) t = 6.3. The right column is for a vortex rotating clockwise 
at (d) t = 2.5, (e) t = 3.8, (f) t = 6.3. The coordinates are 
in units of Lq. The dashed circles mark the position of the 
circular corrals. 



of the vortex. 

The evolution of the vortex in the corrals is shown in 
Fig. 4. The evolution is obtained numerically from the 
GP equation ([1]). After switching off the harmonic trap, 
and with the expansion of the condensate, the condensate 
is reflected by the corral-like potential. Similarly to the 
preceding results, the expanding and reflected conden- 
sates will overlap, and lead to clear interference patterns. 
These interference patterns are clearly seen in Figs. 4a- 
4f, where Figs. 4a-4c show the evolution for the case of 
e*^ while Figs. 4d-4f give the evolution for the case of 

In addition to the sunflower- like structure because of 
discrete characteristic of the corral-like potential, one 
more interesting feature has emerged due to the pres- 
ence of the vortex. At t = 2.5To, as shown in Fig. 4a 
and Fig. 4d, the two different vortices have almost iden- 
tical density distribution even though their sunflower- like 
structure indicates that the two BECs have already felt 
the corral-potential. The rotational difference between 
the two BECs shows up only when the "main peak" of 
the condensate hits the corral-potential, as shown in Figs. 
4b-c and Figs. 4e-f. This is due to the unique property of 
the velocity field for a quantized vortex, where the center 
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FIG. 5: Density distributions of a condensate in a stadium- 
shaped potential at different times: (a) t = 0, (b) t = 0.5T, 
(c) t = T, (d) t = L5T, (e) t = 3T, (f) t = 6T. T = 2{R+d)/v 
is the time needed for a classical particle to propagate along 
the stadium axis. The coordinates are in units of Lq. In Fig. 
(a), the stadium-shaped confinement potential is also shown. 

of the vortex has larger velocity. In Figs. 4c and 4f, a se- 
ries of flyer-like peaks distributed along the circumference 
of the outmost circle distinguish obviously two vortices 
with different rotational directions. This feature may be 
applied to detect the rotational direction of a vortex. 

It is still an open problem to experimentally inves- 
tigate the behavior of a quantized vortex in electronic 
quantum corrals. For an atomic condensate, all the 
rich dynamics of a vortex shown above can be readily 
observed experimentally as a macroscopically quantized 
vortex can be generated experimentally in a BEG with 
mature technologies [l2| 



III. QUANTUM CHAOS FOR A BEG IN 
GORRAL-LIKE POTENTIAL 

We now turn to consider the dynamics of a BEG con- 
fined in stadium-shaped quantum corrals as shown in Fig. 
5a. Stadium billiard is a typical system to study both 
classical and quantum chaos The energy- level distri- 
bution was proposed in Ref. [l3| to reveal the quan- 
tum chaos for a BEG. Most recently, the scars in the 
steady-state density profiles of parametrically driven con- 
densates was proposed to study the quantum chaos pl| . 
Our work here focuses on the dynamic manifestation of 
quantum chaos with a BEG. 

The stadium potential can be described by Voc with 
{xi^yi} uniformly distributing along the circumference 
of the stadium as shown in Fig. 5a. In our numeri- 
cal computation, the parameters that we choose for the 
stadium-shaped quantum corrals are M = 50, a = 1, 
7 = 500, R = 10, and d = 307r/19. The groundstate 
wave function is calculated for the condensate confined 
in the harmonic trap with uj± = 87. bn and the stadium- 



FIG. 6: Density distributions of a condensate in a circle- 
shaped potential at different times: (a) t = 0, (b) t = 0.5T, 
(c) t = T, (d) t = 1.5T, (e) t = 3T, (f) t = 6T. Here 
T = 2R/v. The coordinates are in units of Lq. In Fig. (a), 
the circle-shaped confinement potential is also shown. 



shaped quantum corrals. This condensate is then given 
an initial velocity of v. After switching off the harmonic 
trap, we study numerically the evolution of this moving 
condensate in the stadium-shaped quantum corrals. For 
V = 2.7 mm/s. Figs. 5a-5f show a series of snapshots 
for the condensate evolution. The regular density distri- 
bution is destroyed after a small number of reflections. 
This is in stark contrast with the dynamics of a mov- 
ing BEG in circle-shaped quantum corrals, which stays 
regular even after long-time evolution as shown in Fig. 
6. 

This quantum chaotic behavior was firstly revealed in 
Ref. [8] for a free particle in a continuous stadium bil- 
liard. In this work, Tomsovic and Heller showed that 
the chaotic wave dynamics shown in Fig. 5 can be re- 
garded as the supposition of millions of classical trajec- 
tories. The random looking in Fig. 5 is just a reflection 
of the chaotic classical trajectories. In contrast to a free 
particle, there is interatomic interaction in a BEG and 
the corral potential is discrete. As a result, one might ex- 
pect with slight hesitation that a BEG in stadium-shaped 
quantum corrals exhibits typical quantum chaotic wave 
dynamics. Our numerical calculations above show that 
the hesitation can be cast aside. In fact, no obvious dif- 
ference is found, compared to the situation of noninter- 
acting BEG. We have also considered the situation with 
increasing coupling constant. It is found that increas- 
ing the coupling constant has the effect of enhancing the 
dynamical quantum chaotic behavior. 

As mentioned at the beginning, it has remained elu- 
sive to visualize experimentally the quantum chaotic be- 
havior in the electronic quantum corrals despite a great 
deal of effort [7]. The main reason is that the quantum 
corrals are too "leaky" for electrons. Our above nu- 
merical studies show that the BEG system provides a 
very good chance to experimentally visualize the quan- 
tum chaos: The stadium-shaped corrals can be built with 
the scheme shown in Fig. 1. The BEG can gain a ve- 
locity with the method used for the experimental stud- 
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ies of quantum reflection [iS^. Considering the length 
unit is 5 jam in Figs. 2-6, a imaging resolution below 
5 /im is necessary to reveal the structure in these flg- 
ures. Fortunately, high resolution absorption imaging 
below 4 jam has been achieved ^ISJ The chaotic behav- 
ior has been studied experimentally with nondegener- 
ate ultracold atomic gases [l7j; quantum chaos has been 
visualized with classical wave systems, such as sound 
and miser owaveplsj. An eventual implementation of our 
scheme will be an experimental realization of quantum 
chaos with coherent matter wave, marking an important 
step forward in the experimental studies of chaos. 

Besides visualizing the known quantum chaotic behav- 
ior, we can try to explore more on quantum chaos with 
BECs in a stadium-shaped corrals. For example, we can 
study the thermal effect on quantum chaos by placing 
non-condensed cold atoms and partially condensed cold 
atoms into the corrals. This is impossible for the elec- 
tronic quantum corrals. One can also study how the 
interference between two BECs is affected by quantum 
chaos. We have shown in Fig. 7 the evolution of two 
initially coherently separated condensates [l^. The sig- 
niflcance of this study lies in that a classical particle can 
not be regarded as two coherently separated particles. 
Figs. 7(al)-7(a6) give the evolution of two coherently 
separated Gaussian wavepackets with the same condition 
in Fig. 5. The initial width and distance are ^J^Jh and 
4. For short-time evolution, the interference between two 
condensates is seen in Fig. 7(a2). After long-time evo- 
lution, the interference is destroyed by quantum chaos 
as seen in Figs. 7(a4)-7(a6). As a comparison, we have 
also computed the evolution of these two BECs in circle- 
shaped corrals. The results are shown in Fig. 7(bl-b6), 
where the interference between two condensates is pre- 
served partly even for long-time evolution. In the inset 
of Fig. 7(b5), the fork- like structure implies the informa- 
tion of two initially coherently separated condensates. 
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IV. SUMMARY AND DISCUSSION 

In summary, we have studied an atomic analogy to 
electronic quantum corrals. A scheme is proposed to 
study the dynamic evolution of a Bose-Einstein conden- 
sate conflned in a corral-like potential. In particular, with 
these atomic corrals, it is now promising to study ex- 
perimentally the dynamic quantum chaotic behavior Q. 
The atomic quantum corrals proposed here can also be 
applied to study other quantum gases, such as molecular 
Bose-Einstein condensates [20|, degenerate Fermi gases 
[21], ultracold Fermi gases in the unitarity limit |22], and 
BCS (Bardeen-Cooper-Schrieffer) superfluid [23] in the 
corral-like potential. The different descriptions of the 
order parameter, quantum statistics and evolution equa- 
tion [24] mean that there should be very rich phenomena 
for discovery. 



FIG. 7: The left and right columns give the evolution of 
two initially coherently separated condensates in a stadium- 
shaped and circle-shaped potentials, respectively. Each figure 
in the left and right columns is 30 x 40 and 30 x 30 with co- 
ordinate unit Lq. Here the time unit — 2{R^ d) jv and 
T — 2R/v for the left and right columns. 



Acknowledgments 



We acknowledge useful discussions with Baolong Lii, 
Cheng Chin, and Jin Wang. This work was supported by 
NSFC under Grant Nos. 10875165, 10825417, 10634060, 
and NBRPC 2006CB921406. 



6 



[1] M. F. Crommie, C. P. Lutz, and D. M. Eigler, Nature 
(London) 363, 524 (1993); M. F. Crommie, C. P. Lutz, 
and D. M. Eigler, Science 262, 218 (1993). 

[2] G. A. Fiete and E. J. Heller, Rev. Mod. Phys. 75, 933 
(2003). 

[3] G. C. des Francs et al, Phys. Rev. Lett. 86, 4950 (2001). 

[4] C. Chicanne et a/., Phys. Rev. Lett. 88, 097402 (2002). 

[5] Nature Insight: Ultracold Matter [Nature (London) 416, 
205 (2002)]; J. P. Yin, Phys. Rep. 430, 1 (2006). 

[6] F. Dalfovo, S. Giorgini, L. P. Pitaevskii, and S. Stringari, 
Rev. Mod. Phys. 71, 463 (1999); A. J. Leggett, Rev. 
Mod. Phys. 73, 307 (2001); C. J. Pethick and H. Smith, 
Bose-Einstein Condensation in Dilute Gases (Cambridge 
University, Cambridge, 2002); L. P. Pitaevskii and S. 
Stringari, Bose-Einstein condensation (Clarendon, Ox- 
ford, 2003). 

[7] E. J. Heller, M. F. Crommie, C. P. Lutz, and D. M. 

Eigler, Nature 369, 464 (1994); M. F. Crommie, C. P. 

Lutz, D. M. Eigler, and E. J. Heller, Physica D 83, 98 

(1995); M. F. Crommie, C. P. Lutz, D. M. Eigler, and E. 

J. Heller, Surf. Sci. 361, 864 (1996). 
[8] S. Tomsovic and E. J. Heller, Phys. Rev. Lett. 67, 664 

(1991); Phys. Rev. E 47, 282 (1993). 
[9] S. Tung, V. Schweikhard and E. A. Cornell, Phys. Rev. 

Lett. 97, 240402 (2006). 
[10] D. Clement et a/.. New J. of Phys. 8, 165 (200 6). 
[11] K. Henderson et a/.. Preprint arXiv:0902.2171 ' (2009) ; A. 

Itah et aL, Preprint arXiv:0903.3282 (2009); V. Boyer et 

al, Phys. Rev. A 73, 031402(R) (2006). 
[12] A. L. Fetter, Preprint arXiv:0801.2952 (2008); M. R. 

Matthews et a/., Phys. Rev. Lett. 83, 2498 (1999); K. 



W. Madison, F. Chevy, W. Wohlleben, and J. Dalibard, 
Phys. Rev. Lett. 84, 806 (2000); J. R. Abo-Shaeer, C. 
Raman, J. M. Vogels, and W. Ketterle, Science 292, 476 
(2001). 

[13] C. W. Zhang, J. Liu, M. G. Raizen, and Q. Niu, Phys. 
Rev. Lett. 93, 074101 (2004). 

[14] N. Katz and O. Agam, Preprint arXiv: 0903.0968 (2009). 

[15] T. A. Pasquini et a/., Phys. Rev. L ett. 93, 223201 (2004). 

[16] N. Gemelke et a/.. Preprint 0904. 1532 (2009); pri- 
vate communication with Dr. Cheng Chin about the im- 
proved imaging resolution of about 1.5 /im. 

[17] V. Milner, J. L. Hanssen, W. Campbell, and M. G. 
Raizen, Phys. Rev. Lett. 86, 1514 (2001); D. A. Steck, W. 
H. Oskay, and M. G. Raizen, Science. 293, 274 (2001); 
N. Friedman et a/., Phys. Rev. Lett. 86, 1518 (2001). 

[18] H. J. Stockmann, Quantum Chaos: an Introduction 
(Cambridge, Cambridge, 1999). 

[19] M. R. Andrews et a/.. Science 275, 637 (1997). 

[20] S. Jochim et a/.. Science 302, 2101 (2003); M. Greiner, 
C. A. Regal, and D. S. Jin, Nature 426, 537 (2003); M. 
W. Zwierlein et a/., Phys. Rev. Lett. 91, 250401 (2003). 

[21] B. D. Marco, and D. S. Jin, Science 285, 1703 (1999). 

[22] K. M. O'Hara et a/.. Science 298, 2179 (2002). 

[23] C. A. Regal, M. Greiner, and D. S. Jin, Phys. Rev. Lett. 
92, 040403 (2004); M. W. Zwierlein et a/., Phys. Rev. 
Lett. 92, 120403 (2004); M. Bartenstein et a/., Phys. Rev. 
Lett. 92, 120401 (2004); J. Kinast et a/., Phys. Rev. Lett. 
92, 150402 (2004). 

[24] S. Giorgini, L. P. Pitaevskii, and S. Stringari, Rev. Mod. 
Phys. 80, 1215 (2008). 



